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Abstract. In this paper we discuss variational constrained mechanics (vako- 
nomic mechanics) on Lie aflgebroids. We obtain the dynamical equations and 
the aff-Poisson bracket associated with a vakonomic system on a Lie aflgebroid 
A. We devote special attention to the particular case when the nonholonomic 
constraints are given by an affine subbundle of A. and we discuss the varia- 
tional character of the theory. Finally, we apply the results obtained to several 
examples. 



1. Introduction 

Lie algebroids have deserved a lot of interest in recent years (from a theoretical 
and applied point of view). In the context of Mechanics, an ambitious program 
was proposed by Weinstein [34] in order to develop geometric formulations of the 
dynamical behavior of Lagrangian and Hamiltonian systems on Lie algebroids. In 
the last years, this program has been actively developed by many authors. In fact, a 
Klein's formalism for unconstrained Lagrangian systems on Lie algebroids has been 
discussed and a symplectic formulation of Hamiltonian mechanics on these objects 
has been developed (see [50J [221 I2H 0S3])- The main notion is that of prolongation 
of a Lie algebroid over a map introduced by Higgins and Mackenzie .14. An 
alternative approach, using the linear Poisson structure on the dual bundle of a Lie 
algebroid, was discussed in [T2] , 

An interesting kind of mechanical systems are those subject to external linear 
constraints. For these systems, one may derive the dynamical equations using the 
Lagrange-D'Alcmbcrt principle (nonholonomic mechanics) or using a constrained 
variational principle (vakonomic mechanics). The resultant equations are, in gene- 
ral, different. Constrained Lagrangian systems (variational or not) have application 
in many different areas: engineering, optimal control theory, mathematical econom- 
ics, sub-Riemannian geometry, motion of microorganisms, etc. For a geometrical 
treatment of standard mechanical systems subject to external linear constraints we 
remit to the monographs E] and references therein. 

More recently, several authors discuss the more general class of nonholonomic La- 
grangian (Hamiltonian) systems subject to linear constraints on Lie algebroids (see 
[HI [23 HE] ) ■ In the same Lie algebroid setting, other authors [16] consider varia- 
tional constrained mechanical systems. In another direction, a unified approach of 
nonholonomic and vakonomic mechanics, using general algebroids instead of just 
Lie algebroids, was developed in [5]. 

As a consequence of all these investigations, one deduces that there are several 
reasons for discussing unconstrained (constrained) Mechanics on Lie algebroids: 
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i) The inclusive nature of the Lie algebroid framework. In fact, under the same 
umbrella, one can consider standard unconstrained (constrained) mechanical sys- 
tems, (nonholonomic and vakonomic) Lagrangian systems on Lie algebras, uncon- 
strained (constrained) systems evolving on semidirect products or (nonholonomic 
and vakonomic) Lagrangian systems with symmetries. 

ii) The reduction of a (nonholonomic or vakonomic) mechanical system on a Lie 
algebroid is a (nonholonomic or vakonomic) mechanical system on a Lie algebroid. 
However, the reduction of an standard unconstrained (constrained) system on the 
tangent (cotangent) bundle of the configuration manifold is not, in general, an 
standard unconstrained (constrained) system. 

iii) The theory of Lie algebroids gives a natural interpretation of the use of quasi- 
coordinates (velocities) in Mechanics (particularly in nonholonomic and vakonomic 
mechanics). 

On the other hand, in (TUl US] an affine version of the notion of a Lie algebroid 
structure was introduced. The resultant geometric object is called a Lie affgebroid 
structure. A Lie affgebroid structure on an affine bundle A is equivalent to a Lie 
algebroid structure on the bidual bundle to A such that the section of the affine 
dual to A induced by the constant map on A equal to 1 is a 1-cocycle. 

Lie affgebroid structures may be used to develop a time-dependent version of 
unconstrained Lagrangian and Hamiltonian equations on Lie algebroids (see [111 
[T7J [251 HH1 GH] ) ■ In addition, in [IB] the authors present a geometric description of 
Lagrangian and Hamiltonian systems on Lie affgebroids subject to affine nonholo- 
nomic constraints. If we apply this general theory to the particular case when the 
Lie affgebroid is the 1-jet bundle of local sections of a fibration r : Q — > R then 
one recovers some results obtained in [H 1181 119j for standard time-dependent non- 
holonomic Lagrangian systems subject to affine constraints. The same reasons for 
discussing unconstrained (constrained) mechanics on Lie algebroids are also valid 
for discussing unconstrained (constrained) mechanics on Lie affgebroids. 

On the other hand, in [33) the authors discuss standard time-dependent vako- 
nomic dynamics and its relation with presymplectic geometry. More recently, in 
[I] a geometric approach to time-dependent optimal control problems is proposed. 
This formulation is based on the Skinner and Rusk formalism for Lagrangian and 
Hamiltonian systems. Some applications are also presented. The aim of this paper 
is to extend these formulations to the Lie affgebroid setting or, in other words, to 
discuss vakonomic mechanics on Lie affgebroids. 

The paper is organized as follows. In Section 2, we recall some well-known facts 
about the geometry of Lie affgebroids and about the unconstrained Hamiltonian 
formalism on Lie affgebroids. In Section 3, we obtain the vakonomic equations and 
the vakonomic bracket for a constrained mechanical system on a Lie affgebroid A. 
We devote special attention to the particular case when the constraints are given by 
an affine subbundle of A. We also discuss, in this section, the variational character 
of the theory. In section 4, we apply the results obtained in the paper to several 
examples. In fact, we develop a Skinner-Rusk formalism on Lie affgebroids. We 
also consider vakonomic Mechanics on a Lie affgebroid A, for the particular case 
when A is the 1-jet bundle of local sections of a fibration over E. As a consequence, 
we recover some previous results in the literature. We also discuss optimal con- 
trol systems as vakonomic systems on Lie affgebroids. The paper ends with our 
conclusions and a description of future research directions. 
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2. Hamiltonian formalism on Lie affgebroids 

2.1. Lie affgebroids. Let ta ■ A — > Q be an affine bundle with associated vector 
bundle ry : V — > Q. Denote by r iA + : A + = Aff(A, R) — > Q the dual bundle whose 
fibre over x £ Q consists of affine functions on the fibre A^. Note that this bundle 
has a distinguished section 1a G r(r^+) corresponding to the constant function 1 
on A. We also consider the bidual bundle tj : A — » Q whose fibre at x G Q is the 
vector space A x = (A+)*. Then, A may be identified with an affine subbundle of 
A via the inclusion %a : A — > A given by vi(a)(<£>) = </?(a), which is an injective 
affine map whose associated linear map is denoted by iy : V — > A. Thus, V may 
be identified with a vector subbundle of A. 

A Lie affgebroid structure on .A consists of a Lie algebra structure [•, -]y on the 
space r(ry) of the sections of ry : V — > Q, a R- linear action D : T(ta) x L(ry) —> 
r(ry) of the sections of A on r(ry) and an affine map pj\ : A — * TQ, the anchor 
map, satisfying the following conditions: 

. D X \Y,Z\ V = ID X Y,Z\ V + {Y,D X Z\ V , 

• D x+ yZ = D x Z+{Y,Z] v , 

• D x (fY) = fD x Y + p A (X)(f)Y, 

for X e T(t a ), Y,Z 6 T(t v ) and / e C°°(Q) (see HDHSJ). 

If ([*> *]V) D, pj{) is a Lie affgebroid structure on an affine bundle A then (V, 
IvIvjPv) is a Lie algebroid, where pv ■ V —* TQ is the vector bundle map 
associated with the affine morphism pA ■ A — * TQ (for the definition and properties 
of Lie algcbroids we remit to [22]). 

A Lie affgebroid structure on an affine bundle ta '■ A — * Q induces a Lie algebroid 
structure ([•, -Ja'Pa) on * ne bidual bundle A such that 1a € T(ja+) is a 1-cocycle 
in the corresponding Lie algebroid cohomology, that is, d A lA = 0. Here, d A is the 
differential of the Lie algebroid (A, [•, -Ix-iPa)' 

Conversely, let (U, [•, -JutPu) be a Lie algebroid over Q and <j) : U — > K be a 
1-cocycle of (17, [•, pj/) such that 0^ ^ 0, for all x € Q. Then, A = _1 {1} 
is an affine bundle over Q which admits a Lie affgebroid structure in such a way 
that (U, [•, -\u,Pu) ma Y be identified with the bidual Lie algebroid (A, [•, 'I?, Pa) 
to A and, under this identification, the 1-cocycle 1a '■ A — > R is just </>. The affine 
bundle t/l : A — > Q is modelled on the vector bundle ry : V — (/> _1 {0} — > Q. 

Let T4 : A — > Q be a Lie affgebroid modelled on the Lie algebroid ry : V — » Q. 
Suppose that (a; 1 ) are local coordinates on an open subset {/ of Q and that {eo, e a } 
is a local basis of r(rr) in J7 which is adapted to the 1-cocycle 1»., i.e., such that 
1-4 (eo) = 1 and l^ec) = 0, for all a. Note that if {e°,e a } is the dual basis of 
{eo, e Q } then e° = 1a- Moreover, since 1a is a 1-cocycle, we have that 

■ d d 
[e ,e Q ]^ = C y Q; e 7 , [e Q , = p^(e ) = p ^— , P^(e a ) = p* a ^ . 

Denote by (x l ,y a ,y a ) the corresponding local coordinates on A. Then, the local 
equation defining the affine subbundle A (respectively, the vector subbundle V) 
of A is y° — 1 (respectively, y° = 0). Thus, (x l ,y a ) may be considered as local 
coordinates on A and V. 

The standard example of a Lie affgebroid is the 1-jet bundle ti.o : J 1 t — > Q of 
local sections of a fibration r : Q — > R. It is well known that ti : o is an affine bundle 
modelled on the vector bundle 7r = (7TQ)|y r : Vt — > Q, where V^r is the vertical 
bundle of r. Moreover, if t is the usual coordinate on R and is the closed 1-form 
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on Q given by r\ — T*{dt) then we have the identification J t = {»£ TQ \ r)(v) = 1} 
(see, for instance, [32). Note that Vt — {v G TQ\rj(v) = 0}. Thus, the bidual 
bundle J 1 t to T\ o : J t — > Q may be identified with the tangent bundle TQ to Q 
and, under this identification, the Lie algebroid structure on ttq : TQ — > Q is the 
standard Lie algebroid structure and the f-cocycle lji T on 7Tq : TQ — > Q is just rj. 

2.2. The Hamiltonian formalism. Suppose that (tji : A — > Q, ry : V — > Q) is a 

Lie afFgebroid. Then, we consider the prolongation T 71 ^* of the bidual Lie algebroid 
(A, [•, -]^, p^) over the fibration t v : V* — ► Q and denote by ([•, ■Jjr jPjr ) tne Lie 

algebroid structure on T^y* (for the definition of the Lie algebroid structure on 
the prolongation of a Lie algebroid over a fibration, we remit to [TH |2T)] ). 

Let \i : A + — > V* be the canonical projection given by fi(ip) — <p l , for ip G .A+, 
with x € Q, where (/r G V* is the linear map associated with the affine map ip and 
h : V* — > .A + be a Hamiltonian section of /i, that is, /i° h = Id. 

Now, we consider the prolongation T^yi" 1 " of the Lie algebroid j4 over r A + : j4 + — > 
Q with vector bundle projection : 7 A A + — > ,A + . Then, we may introduce the 

map T7i : T^l/* -> 7 A A+ defined by 07i(a,X a ) = (&,(T a h)(X a )), for (a,X Q ) G 
T^V*, with a E V* . It is easy to prove that the pair (7h,h) is a Lie algebroid 
morphism between the Lie algebroids r-' : T^y* — > V* and t^ 4-1- : 7 A A + — > .A + . 
We denote by A^ and Q/j the sections of the vector bundles (7 A V*)* — » y* and 
A 2 (T-*V*)* -> y* defined by 

A ft = (T/i,ft)*(A^), n^^/i)*^), (2.1) 
A j- and f2 jj being the Liouville section and the canonical symplectic section, respec- 
tively, associated with the Lie algebroid A (see (5U]). Note that flh = — d 7 v A/,. 
On the other hand, let 77 : 7 A V* — ► M be the section of (T^V*)* — > y* given by 
77(5,^) = U(a), for (a,X v ) G 7 A V* , with 1/ G V*. (2.2) 

We remark that if pi\ : 7 A V* — > yi is the canonical projection on the first factor 
then (pr 1; Ty) is a morphism between the Lie algebroids t-' : 7 A V* — > V* and 
: j4 — > Q and (pr 1; r{>)*(lyi) = 77. Thus, since 1^ is a 1-cocycle of : A — » Q, 
we deduce that 77 is a 1-cocycle of the Lie algebroid r~ v : T^y* — > V*. 

Let (a; 2 ) be local coordinates on an open subset U of Q and {eo,e a } be a local 
basis of r(rr) on /7 adapted to 1a- Denote by (x l , y a , y a ) the induced local coordi- 
nates on A and by (a; 1 , yo,y a ) the dual coordinates on A + . Then, (x l , y a ) are local 
coordinates on V* and {^0, V Q , U a } is a local basis of r(r~ / ), where 

y„M = (eo(xU^ J^W = (e^pi^jM) = (0, A J, 

for V G V£. Suppose that h(x\y a ) = (x\ -H(x j ,yp),y a ) and that {y ,^,!! } is 
the dual basis of {^0,^0, W Q }. Then 7/ = y° and, from (|2.ip and the definition of 
the map 7h, it follows that 

n h = TAU a + -cr^T a ^ + (p^— - c2 aVl )T a y° + — u a a y°. 

Thus, it is easy to prove that the pair (Q^jV) i s a cosymplectic structure on the 
Lie algebroid t t $ : 7 A V* y* (this means that d^^fl,, = 0, d 7 "^*/? = and 



CONSTRAINED MECHANICS ON LIE AFFGEBROIDS 



5 



(rjAfth A. .." ■■■An h )(ij}) ^ 0, £ V*, where n is the rank of A) . If R h £ r(r~ v ) 
is the Reeb section of (Qh,v) (that is, in h £lh = and ij^Jj = 1), then its integral 
curves (i.e., the integral curves of p T ~{Rh)) are just the solutions of the Hamilton 
equations for h, 

dx l _ i i dH dy a _ i dH / 7 7 <9i/\ 

for i £ {1, . . . , m} and a £ {1, . . . , n}. 

Next, we will present an alternative approach in order to obtain the Hamilton 
equations. For this purpose, we will use the notion of an aff-Poisson structure on 
an AV-bundle which was introduced in [10] (see also [11])- 

Let tz : Z — > Q be an affine bundle of rank 1 modelled on the trivial vector 
bundle tq x r :QxR^Q, that is, tz : Z — > Q is an AV-bundle in the terminology 
of [11] . Then, we have an action of R on the fibres of Z. This action induces a 
vector field Xz on Z which is vertical with respect to the projection tz '■ Z — > Q. 

On the other hand, there exists a one-to-one correspondence between the space 
of sections of r z : Z -> Q, T(r z ), and the set {F h £ C°°(Z) | X z (F h ) = -1}. In 

d 

fact, if /i £ r(rz) and {x x , s) are local fibred coordinates on Z such that Xz = tt- 

ds 

and h is locally defined by h{x l ) = (x l , — H (x 1 )) , then the function i 7 ^ on Z is 
locally given by Fh{x l , s) = —H(x % ) — s, (for more details, see [llj). 

Now, an aff-Poisson structure on the AV-bundle tz ■ Z — > Q is a bi-affine map, 
{■, ■} : r(rz) x r(rz) — > C°°(Q), which satisfies the following properties: 

i) Skew-symmetric: {hi, /12} = — {h,2, h\}. 

ii) Jacobi identity: {hi, {h 2 , h 3 }} v + {h 2 , {ha, h{\} v + {ha, {hi, h 2 }}v = 0, 
where {•, -}y is the afhnc-lincar part of the bi-affine bracket. 

iii) If he T(t z ) then the map {h, ■} : T{t z ) -> C°°(Q) defined by {h, -}{h') = 
{h,h'\, for h! £ T(tz), is an affine derivation. 

Condition iii) implies that, for each h £ T(rz) the linear part {h, -}y : C°°(Q) — > 
C°°(Q) of the affine map {ft, •} : r(rz) — > C 00 ((5) defines a vector field on Q, which 
is called the Hamiltonian vector field of h (see 

In [TT], the authors proved that there is a one-to-one correspondence between 
aff-Poisson brackets {•, •} on tz : Z — > Q and Poisson brackets {•, -}n on Z which 
are X^-invariant, i.e., which are associated with Poisson 2-vectors II on Z such 
that £jx z II = 0. This correspondence is determined by 

{hi,h 2 }°T Z = {F hl ,F h2 } n , for hi,h 2 £ T(tz). 

Using this correspondence one may prove the following result. 

Theorem 2.1. [17j Let tj\ : A — > Q be a Lie affgebroid modelled on the vector 
bundle Ty : V — > Q. Denote by tj±+ : A + — > Q (resp., Ty : V* — > Q) the dual vector 
bundle to A (resp., to V) and by /j : A + — > V* the canonical projection. Then: 

i) (j, : A + — > V* is an AV-bundle which admits an aff-Poisson structure. 

ii) If h : V* — > yi + is a Hamiltonian section then the Hamiltonian vector field 
of h with respect to the aff-Poisson structure is a vector field on V* whose 
integral curves are just the solutions of the Hamilton equations for h. 

3. Vakonomic mechanics on Lie affgebroids 

3.1. Vakonomic equations and vakonomic bracket. Let ta ' A — ► Q be a Lie 

affgebroid of rank n over a manifold Q of dimension m. We consider an embedded 
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submanifold M C A, called the constraint submanifold, of dimension n + m — m 
such that -rjvt = t a \ m : M — > Q is a surjective submersion. 

Now, suppose that e is a point of M, with 7>t(e) = x, that (a; 2 ) are local coordi- 
nates on an open subset U of Q, x eU, and that {eo, e a } is a local basis of T(tj) 
on 17 adapted to the 1-cocyclc Denote by (x l ,y°, y a ) (respectively, (x l ,y a )) the 
corresponding local coordinates for A (respectively, A) on the open subset T j^{U) 
(respectively, r^ 1 ([/)). Assume that 

JVC n t~\U) = {(x\y a ) G t a \U) | ^ A (x\y a ) = 0, A = 1, . . . , m}. 

The rank of the (m x (n + m))-matrix ( — — -, — — is maximun, that is, m. Then, 
v v ;; V dx l dy a ) 

using that tj^ : M — > Q is a submersion, we can suppose that the (m x m)-matrix 

V dy B e / A=l,...,m;S=l,...,m 

is regular. Then, we will use the following notation (y a ) = (y A ,y a ), for 1 < a < n, 
1 < A < rh and rh + 1 < a < n. 

Now, using the implicit function theorem, we obtain that there exist an open 
subset V of t^(U), an open subset W C R m+ ™ m and smooth real functions 
* A : W -> K, A = 1, . . . , m, such that 

Mn^E {(a;*, y") e V- 1 y A = * A (x\y a ), A = l,...,fh}. 

Consequently, (x l ,y a ) are local coordinates on M. 

Next, consider the Whitney sum of A + and A, that is, A + ®qA and the canonical 
projections p^ : A + ®q A — > and pr 2 : ,A + ©q ./I — > A Let Wo be the 
submanifold of A + ®q A given by Wo = pr~ (M) = (Bq M and the restrictions 
7Ti = pri|vi/ an( i ^2 = P r 2|wv ^lso denote by v : Wo — > Q the canonical projection. 

Now, we take the prolongation , r~' t+ : 7 A A + — > ,A + (respectively, r~ : T^Wo — > 

Wo) of the Lie algebroid A over t iA + : A + — > Q (respectively, f : Wo — > Q). 
Moreover, we can prolong 7Ti : Wo — ► to a morphism of Lie algebroids Ttti : 
T A W -> 7 A A+ defined by Ttti = (Id, Tin). 

If (x l ,yo, y a ) are the local coordinates on A + induced by the dual basis {e°, e a } 
of the local basis {eo, e a } of T(tj-), then (x l , yo,y a ,y a ) are local coordinates for Wo 
and we may consider the local basis {^o, ^a, 3**, V a } of r(r~) defined by 



A' 

%(ip,a) = (e o (x),p o -^-r^,0y y a (ip,a) = (e a (x), > o) , 



T°(^a)= (O,/- , ), ? Q (^,a)= (0,-^- ,o), V„( V ,a) = (o,0, ^ ), 

for (ip, a) G Wo and f(tp, a) = a;, where p and p^, are the components of the anchor 
map with respect to the local basis {eo, e a }. 

Now, one may consider on the Lie algebroid r~ : 7 a Wq — ► Wo the presymplectic 

2-section f2 = (T7Ti, 7Ti)*0^, where fi^ is the canonical symplectic section on T^-A" 1 ". 
The local expression of is 

n = y° a To + T a y Q + c 7 Q y 7 y° a y« + ^^y a ^, (3.1) 

{^^".To^a.V } being the dual basis of the local basis y a , 9°, V a }. 
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On the other hand, if pr x : 7 Wo — > A is the canonical projection on the first 
factor, then we can introduce the section r\ £ T((t~)*) defined by 77 = (pr 1; v)*1ji. 

Since 1a is a f-cocycle of A — > Q, we deduce that 7/ is a f-cocycle of T^Wo — > Wo- 
Moreover, it is easy to prove that 

r? = y°. (3.2) 

Now, let L : A — > K be a Lagrangian function on 71 and denote by L the 
restriction of L to the constraint submanifold M. 

The Pontryagin Hamiltonian Hw is the real function in Wq — A + ©q M given 
by Hwoi'fi, a) = <^(a) — L(&), or, in local coordinates, 

H Wo (x\yo, y a ,y a ) = yo + y a y a + y A ^> A (x\y a ) - L{x\y a ) . (3.3) 

Thus, one can consider the presymplectic 2-section f2vi/ on 7 a Wq defined by 

n Wo = n + d 7 " 1 w ° H Wo At?. 

In local coordinates, using (|3.ip . (|3.2| and (|3.3| . we deduce that 



fiw„ = T a y a + 



( 



d^ J 



dx 1 dx 



dL\ i 7 



r a y° + a t 



+* A ?A A y° + -C> 7 y° A + (y a + y" - £-) V A 



(3.4) 



Definition 3.1. The vakonomic problem (L,M) on the Lie affgebroid A consists 
of to find the solutions for the equations 

ix^Wo — an d ixV — 1, with X 6 T(t~). (3-5) 

First, we will obtain the local expression of the vakonomic problem. In general, a 
section X satisfying the equations (|3.5p cannot be found in all points of Wo- Thus, 
we consider the points where ([3.5)1 have sense. We define 



W x = {w e W I 3 Z e 7 A W : i Z n Wo (w) = and i zV (w) = 1}. 
In local coordinates, we deduce that W\ is characterized by the equations 



fa = y a + yA' 



d^ J 



dL 



QyO, QyO, 



= 0, m + 1 < a < ri. 



Moreover, a direct computation, using (|3.2|) and ([3.4p , proves that the local 
expression of any section X satisfying the equations ()3.5[) is of the form 

d*f A < 



X 



(T ,T») 



% + * A y A + y a y a + t ? u + 



- yi {Cl +* A Cl A + y*CZ a ) 



Pa 



dL 



\dx l 



y A 



dx 1 



rpa + T a VQ; 



with To and T a arbitrary functions. Consequently, the vakonomic equations are 

Po + V A P i A + y a pi, 



yA 



( dL 
\dx % 

dL 



yB 



dx 
d^ A 



T )pA-y,(cio + ^ B ci B + y a ci a ), 



dt { dy a VA dy a 



( dL 

\dx l 



yA- 



d^ / 

dx 1 



(3.6) 



Pa 



-y,(C2o + ^ B ci B + y b c2 b ), 

for all 1 < i < to, 1 < A < m and to + 1 < a < n. 
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Remark 3.2. The motion equations for the vakonomic mechanics may be also 
expressed as follows 



x=p' 
d 
dt 



y a P l a , 



d_ ( dL_ \ . dL_ _ y 
jAOy") Pa Qx*~ A 



' d /Q<p 

-Jt\d~y a ) ya dx* 



dy a 



(3.7) 



b A = 

.A 



— \$ A an d \ A — y A — J^. Note that in contrast to equations 



where <j) A 

equations (|3.7|) are expressed in terms of the global Lagrangian L : A — > R. Thus, 
the equations Q3.6P stress how the information given by the Lagrangian L outside 
M is irrelevant to obtain the vakonomic equations. This is in contrast with what 
happens in nonholonomic mechanics (see |15|). o 

Then, we know that there exist sections X of T^WolWi — * satisfying (|3 . 5[) . 
However, X doesn't belong, in general, to 7 A W\ C A x TW\. In fact, one may 
prove that the restriction to W\ of Xcy x ») is a section of T^VKi — > W\ if and only 
if 

[T a (<f"^ W °(pii)(V a ) = (d^ W °ip b )(r a V a -X ( y 0! y a) )} ]Wl , V6. 

Then we have a system of (n — m) equations with (n — fn) unknowns (the functions 
T a ). Thus, if we denote by ft a b and (it, the functions 



ft„ 



[(d rAw ^ b )(V a )] lWl 



d 2 L 



yA 



d 2 ^ A 

QyaQyb 



\Wi 



, dy a dy b 

to = [(^ w V 6 )(T«V o -X (T0!Xo) )]| Wl , 

it is clear that the above system has a solution T a if the matrices ft = (ft a b) and 
ft,, = (ft a fc; Mb) have the same rank. Note that even if the above system has a unique 
solution (i.e., if the matrix ft = (ft a f>) is regular), the solution (X(y q ,t°)) i Wi is n °t, 
in general, unique (since the function (To)|wi * s s tiU arbitrary). 

To solve the above problem, we consider a suitable submanifold W[ of W\ whose 
intrinsic definition is 

W[ ={w&W l \ H Wl (w) = 0}, 

where Hw x '■ W\ — > K is the restriction to W\ of the Pontryagin Hamiltonian Hw ■ 
In local coordinates, the submanifold W[ is given by the equation 

yo + y A * V> y h ) + y a y a - U*', y b ) = o. (3.8) 

Let f2w (respectively, ryw^) be the restriction of 17vv (respectively, 77) to 7 A W[. 
Note that the restriction v[ : W[ — > Q of ^ : Wo — > Q to W{ is a fibration and, 
therefore, we can consider the prolongation 7 A W[ of the Lie algebroid A over 
Moreover, we have the following result. 

Proposition 3.3. (fij^'j^W") * s a cosymplectic structure on T A W{ if and only if 
for any system of coordinates (x l , yo, y a , y a ) on Wq we have that 



det(ft afc ) = det 



2 L 

QyaQyb 



d 2 ^ A \ 

VA dy aQyb j ^ °j f° r al1 Voint in W[. 



Proof. It is clear that cP^O^ = and dP Aw '^rtw[ = 0. 
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Now, suppose that the matrix (3l a &) is regular. Since the rank of the Lie algebroid 
7 A W[ — ► W[ is (2n + 1), we have to prove that 

kern W /(ioi) nker?7^K) = {0}, Vwi G W/. 
Now, let Z G kerfi^'C^i) ^ ker rjw (w[ ) . From (|3.4p . it follows that 

(z z r!w„K))(5' K)) = (zz^ K))(VaK)) = 0, for all a. 
On the other hand, 

(d T ^°tf Wo )K)CP°K)) = 1, (d T " V Vfc)K)(VaK)) = for all 6. 

Thus, 7°{w[) £ 7 A ,W[, V a (w[) £ 7 A ,W{ and Z G ker^oK) n kerjjK). This 
implies that Z = A 5 )0 (w / 1 ) + A Q V a (>i). Therefore, since Z G T£, W{, we have that 

= A (^ Wo ^)K)(T°K)) + A°(tP* w V&)(K)(V«K)) = A a 3UK), 
for all &, and, consequently, A a = 0, for all a. Thus, Z = AoCP°(w' 1 ) and 

= Xoid^^'HwMKy^w',)) = A , 

that is, Z = 0. 

The converse is proved in a similar way. □ 



Remark 3.4. We remark that the condition det (3i Q b) ^ implies that the matrix 
d<p a \ 

t I is regular. Thus, using the implicit theorem function, we de- 

&y J a,b=m+l,...,n 

duce that there exist open subsets Wq C Wq, W C R m+ ™ and smooth real functions 
/i a : W — > R, a = m + 1, . . . , n, such that W\ (~l Wo is locally defined by the equations 

y a = ii a (x\y a ), a = m + l,...,n. 

Therefore, we may consider {x l ,yo,y a ) as local coordinates on W\ and, conse- 
quently, from (|3.8p . we obtain that (x l ,y a ) are local coordinates on W{. Thus, a 
local basis of sections oi7 A W{ — > W( is given by {Voi'i Vai'i ^i/}) where 

9Z „.. 9 ^W , „^ M a _ v 



9x l / iWj' 

®/ J,a w \ rna _ (rna ..amO , 9/J, b 



V <9?m /iw/ V dy a J\w[ 



Proceeding as in the proof of Proposition 13.31 we deduce the following result. 

Theorem 3.5. If ,rjw) is a cosymplectic structure on the Lie algebroid 
T i . •j-A-ffli^ jy/ ^/j ew ^/j ere exists a unique section £i G r(r~ 1 ) solution of the 
vakonomic problem (L,M). /n /act, (i is t/ie i?ee& section of (Clyy ,r)w>), that is, 
£i is characterized by the conditions i^Clw^ = and iciVwi = !• 

The above results suggest us to introduce the following definition. 

Definition 3.6. The vakonomic system (L,M) on the Lie affgebroid tj\ : A — > Q 
is said io 6e regular i/ ifte pair {^w^iVwi) * s a cosymplectic structure on the Lie 

algebroid t- 1 : 7 A W{ -> W{. 
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In what follows, we will suppose that (L, M) is a regular vakonomic system on 
the Lie affgcbroid A. Then, from Theorem l3.5[ we have that the vakonomic problem 
has a unique solution which is the Reeb section £1 of the cosymplectic structure 

First, we will give the local expression of the solution section £1. Suppose that 
(x l ,y a ) are local coordinates on W[ as in Remark l3~4l and that {^oi' > ^aV , ^y} is 
the corresponding local basis of T(t~). Then, if {^\, , V" , y a v} is the dual basis of 
{yoi'^ai',^}, we have that (see ([3T3])1 



+ M^al'Ay?, 
Thus, we obtain that 



Pa [va 



dx l 



a(x j ,yp) = %v+V a {xJ,yp)y al ,+y A (x3,n a {xi,yp))y A y 

2/7 (clo + ^ A {x^^{x=,y p ))Cl A + ^{x^yp)Cl a ) 
d^ A dL 



(3.9) 



P a {vA 



3x l \(xi,^(xi,y )) dx l \( x i 4 L°-(xi, Vl3 )) 



7y, 



Now, we will introduce an aff-Poisson bracket on the AV-bundle determined by 
the constraint submanifolds W-y and W{. For this propose, we define the application 
/ii : W\ — > W{ given by 

H\ (ip, a) = (ip - H Wl (ip, a) I A (x) , a) , 

for (cp, a) G W\ C Wq — A + ©q M, with v x {ip, a) = x G Q. 

If (x 1 ,yo,y a ) (respectively, (x z ,y a )) are local coordinates on W\ (respectively, 
W{) as in Remark |3.4[ we deduce that the local expression of fix : W\ — > W[ is 

Pi(x l ,y ,y a ) = (x\y a )- 

Moreover, we have the following result. 

Theorem 3.7. // (L,M) is a regular vakonomic system on the Lie affgebroid tji : 
A — > Q, then fii : W± — > Wj is a AV-bundle which admits an aff-Poisson structure. 

Proof. It is easy to prove that \x\ : W\ — > W[ is an AV-bundle (see Section [32]) . In 
fact, if w = (ip, a) <G (Wi)^, with x E Q, and t G R then 

w + 1 = {ip, a) + 1 = (</5 + tlyi(a;), a). 

To define an aff-Poisson bracket on /Ltj we will introduce a Poisson bracket on W\ 
which is invariant with respect to Xyy t ■ Here, X\y t is the infinitesimal generator of 
the principle action of M. on W\ . 

Consider the prolongation 7(ir\)\Wi '■ f Wi —* T A A + of the restriction (tti)^ : 
W\ — > A + to W\ of the application ni — pxuw '■ Wo - * A + . It is clear that 
{7{'iti)\Wi>( 7r i)\Wi) is a Lie algebroid morphism and, therefore, we can introduce 
the 2-section Q Wl G T(A 2 (t^)*) defined by 



fij- being the canonical symplectic 2-section on 7 A A + . Obviously d 7 Wl Qw! = 0. 
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If (x 1 , yo, y a ) are local coordinates on W\ as in Remark [ITU we can consider the 
local basis of sections {Vol, Vai, ??} of 7 A Wi -> Wi given by 

Voi - (yo + 4gv ) |m , y al = („ a + ^v.) im , 
n = ( o>° Vl , tj = (a» + ^v )^. 

If {y?,y?,T i,T Q i} is the dual basis of {yoi^ai, 3>i, J*" }, we obtain that 

n Wl = y? a ?oi + y? a y al + c^y? a y? + \c\y-Xi a yf . 

Thus, we deduce that flwi is a symplectic section of T 71 W\ — > Wi and, therefore, 
it induces a Poisson bracket on Wi which is given by 

{F, G} Wl = Q Wl (IK^ 1 , M% Wl ) - P J (Mg Wi fOT F,G £ C°°(Wt), 

where 5£p Wl and ^K^ Wl are the Hamiltonian sections associated with the functions 
F and G, respectively, with respect to the symplectic structure Slwi ■ 

Moreover, the Poisson 2- vector il^i determinated by the bracket {•, -}wi is in- 

d 

variant with respect to Xyy, ■ In fact, we have that Xw. — n — an d 

dyo 



tt { d d i d d 

n ^ = ^o^A — + Pa — A — 

7 d d 1 7 9 9 
-^2/7^ A — - 2^^- A — . 



(3.10) 



Thus, we conclude that : W\ — > W x admits an aff-Poisson structure which we 
denote by {•, -} U afc : r(^i) x r(Mi) - * C°°{W[). This structure is characterized by 
the following condition 

{ti 1 ,h1} vak ° l i 1 = {F K ,F h ,,} Wl , for ftU;'er( Ml ), 

•^V ■ -Ffc" being the real functions on W\ associated with the sections h[ , h" (as we 
know, jTwiCFfc/) = ^(Ffc") = -!)■ 

□ 

The aff-Poisson bracket on the AV-bundle pi : W\ — ► W{ , 

{vW:rHxr( Ml )^n^), 

is called £/ie vakonomic bracket associated with the regular system (L,M). 

On the other hand, note that the restriction Hv/x to W\ of the Pontryagin 
Hamiltonian Hw verifies that Xyy 1 {—Hw 1 ) = — 1. Therefore, there exists hi 6 
r(/ii) such that i*^ = —Hw\ - In fact, /ii is the inclusion of into Wi. Moreover, 
we have the following result. 

Theorem 3.8. If F{ G C oa (W[) then the temporal evolution of F[, F[, is given by 

F[ = {hi,F{}" ak , 

{•, ■}^ l ak being the afftne-linear part of the bi-affine bracket {•, -} va k- In other words, 
the Hamiltonian vector field associated with hi with respect to the vakonomic bracket 

coincides with p 1 - (£i ) . 
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Proof. We know that the Hamiltonian vector field {hi, ■}^, l ak of hi with respect to 
the vakonomic bracket is given by 

{hi, ° 1*1 ={Fk 1 ,<P°I*}w 1 , for ^r(lf|). (3.11) 

Then, from ()3.9|) . p.lOp . (|3.1ip and Remark [33 we deduce that this vector field is 
just p-HCi) (see 038). □ 

Next, let /i" be two sections of /ii : Wi — » W{ and suppose that 
h[(x l ,y a ) = {x\-H[{x 3 ,y ),y a ) and h'{(x\y a ) = (x\ -H'{ (x j , y p ), y a ). 
Then, using (I3.10|) . we have that 



{hi,h x Uak - Po + M~9^~9^~ d^d?" 

d(H[-H[>) dH[dH'{ 
" 02/7 dy a C ^ dy a dy ' 



(3.12) 



Since A is a Lie affgebroid, it follows that A + is the total space of an AV-bundle 
over V* with projection fi : A + — > and, moreover, the linear Poisson structure 
11^+ on A + (induced by the Lie algebroid structure of A) defines an aff-Poisson 
bracket {•,•} : T((jl) x r(/x) -> C°°(V*) on the AV-bundle fi : A+ -> V* (see 
Theorem O). 



On the other hand, we may consider the applications (^i)^ : W\ — > A + and 
M ° (^l ) | w[ '■ W{ V* and it is clear that [i°(iri)\ Wl — n° (tti)\w' ° 

In fact, using (|3 . 1 0[) and Remark l3.41 we can prove the following result. 

Corollary 3.9. If (L, M) is a regular vakonomic system on A, then the pair 
((711)1^, /i° (tti)\w') * s a local aff-Poisson isomorphism of AV-bundles, that is: 

i) {tti)\w{ '■ W[ — > V* is a local diffeomorphism; 

ii) The restriction of (^1)1^ to each fibre of fii : Wi —> W{ is an affine 
isomorphism over the corresponding fibre of fi : A + — > V* and 

iii) Ifh'ijh" £ r(/ii) and h',h" G T(/i) satisfy that 

{Tri)\ Wl ° h 'l = h' ° H°{-Kl)\ W ^ {^\)\Wi°K = /l "°M°( 7r i)|iv 1 s 

tterc {/ii,/i'i}wafc = { /l ' ! /l "} A 1 °( 7r i)|w 1 '- 

Remark 3.10. If h°(tti)\w[ '■ W[ — > is a global diffeomorphism then we can 
define the section /i £ T(/i) given by 

h = (iri)\ Wl °hi°(iJ,°(7Ti) lw ,)~ 1 

and it is clear that (Tri)^ °h\ — h°fi° {tti)\w[- This implies that the Hamiltonian 
vector fields of hi and h are {{^i)\Wi , \x ° (7Ti)|w)-related. Therefore, if 7^ : I — > W{ 
is a solution of the vakonomic equations for the system (L, M), then fi ° (tti)\w' ° 7i : 
/ — > 1/* is a solution of the Hamilton equations for /i. Conversely, if 7 : / — > V^* 
is a solution of the Hamilton equations for h then (/i° (7Ti)ny') _1 °7 : — > W[ is a 
solution of the vakonomic equations for the system (L, M). o 

Remark 3.11. If (L,M) is a vakonomic system on a Lie affgebroid j4 which is not 
regular then one may apply a constraint algorithm in order to obtain solutions of 
the vakonomic equations in a suitable Lie subalgebroid of 7 A Wi — > W\. o 
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3.2. Mechanical systems subject to affine constraints on Lie affgebroids. 

The Lagrangian function L : A —> R of a mechanical system on a Lie affgebroid 
ta : A — > Q is of the form 

L(a) = ^Srx(a)(u(a),i>i(a)) - V(r>i(a)), for all a e A, 

where S is a bundle metric on A and V a function on Q. We also denote by S the 
bundle metric induced on A + and we suppose that the 1-cocycle Iji has constant 
norm equal to 1. Moreover, using the metric, we can identify A and V. In fact, we 
have the affine bundle morphism "J : A ^ V given by 

3(a) = iyi(a) — l^(r/i(a)), for all agi, 

where 1 A G T(r^) is defined by tp(l A (x)) = Sx{1a{%), v)j for all ^ S .A^ . Then, if 
S is the restriction to V of S, the Lagrangian L may be written as follows 

£(a) = ^Sr^(a)(J(a),a(a)) - V(r^(a)), for all aei, 

where V(x) = V(x) — |, for all x G Q. 

Now, let (x*) be local coordinates on an open subset U of Q. Then, since the 
section lj\ has constant norm equal to 1, we can consider an orthonormal basis 
of sections of the vector bundle t a \(U) — > U of the form {e° = lyi,e a }. Thus, 
its dual basis {eo, e a } is an orthonormal local basis of sections of A. Moreover, if 
(x l , y°,y a ) are the corresponding local coordinates on A, the local expression of the 
Lagrangian function is 

L{x i ,y a ) = \{y a ) 2 -V{x i ) 

and the Euler-Lagrange equations (that is, the vakonomic equations for the system 
(L,A)) reduce to 

dt -p» + p «y > dt ~ Pa dx> ^o + c a fi v )v ■ 

Next, suppose that the constraint submanifold M of the vakonomic system is an 
affine subbundle 23 of A, that is, we have an affine bundle 23 over Q with associated 
vector bundle tjj% '■ U% — > Q and the corresponding inclusions is : 23 — > A and : 
U'R — > V. Furthermore, assume that \\ G T(t b ). Then, we can choose an special 
coordinate system adapted to the structure of the problem as follows. In fact, 
we consider local coordinates (x l ) on an open subset U of Q and an orthonormal 
local basis of r(ry), {eA,e a }, adapted to the decomposition V = © U%, 

being the orthogonal subbundle to U-b with respect to the bundle metric $. 
Thus, we deduce that {1^ = eo,eA,e a } is an orthonormal local basis of r(r^) 
adapted to the affine subbundle 23. Denote by (x 4 , y°, y A 1 y a ) the corresponding 
local coordinates on A and by (x 1 ,yo,yA,ya) the dual local coordinates on A + . 
Note that the equations which define to 23 as an affine subbundle of A are y A = 0. 
Therefore, the vakonomic system (L, 23) is regular and the local expression of the 
vakonomic equations is 



dt 

dy a _ _ i dV_ 
dt ~ Pa dx l 

y° = l,y A = 0, y a = y a , y = \{y a f + V(x 4 ). 



~Pa~xU ~ (^lo + C2. a y a )y~f, 



2 
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3.3. The variational point of view. Let tj\ : A — > Q be a Lie affgebroid modelled 
on the Lie algebroid ry : V — > Q and L : A — > M be a Lagrangian function on A. 
Next, we will show how to obtain the Euler-Lagrange equations on the Lie affgebroid 
A from a variational point of view. 

We define the set of .A-paths as follows 

Adm([t ,ti],A) = {a: [t ,*i] -^A\p A °& = ^(t a °&)}, 

that is, as the set of admissible curves in A. Then, for two fixed points x,y G Q, 
denote by Adm([t , ti],A)% the set of A-paths with fixed base endpoints equal to x 
and y. 

Now, if iy : V — > A is the canonical inclusion, we consider as infinitesimal 
variations the complete lifts of sections of Ty : V — > Q which vanish at the points 
x and y, that is, 

{(iv°X)\ A I * e r ( r v), X(ar) = and X(y) = 0}. 

Note that if {e ,e a } is a local basis of L(r^) and X G r(ry) is locally given by 
X = X Q e a , then (iy °X)? A is the vector field on A given by 

where = X>^, X- = ^(p* + /p^) - X^(C« + /C^), for all i and a. 

On the other hand, we introduce the action functional SS : Adm([t 0} ti],A) — > M 
defined by 

55(a) = / 1 L(&{t))dt. 

J t 

With this definition it is not difficult to prove that the critical points of SS on 
Adm([to,ti],A)% are the curves a <G Adm([to,ti],A)% which satisfy the Euler- 
Lagrange equations (that is, the vakonomic equations for the system (L,M), with 
M = A). 

Now, let (L,M) be a vakonomic system on the Lie affgebroid t a : A — > Q. 
Denote by Ad77j([t , ii], M)| the set of A-paths on M with fixed base endpoints 
equal to x and y, respectively, that is, 

Adm{[t ,ti\M)l = {& e Adm{{t Q ,t 1 ],A) y x \&(t) eM,Vte [to,h]}- 

In this case, we are going to consider infinitesimal variations (that is, complete 
lifts (iy a X)° A , with X G L(tv)) tangent to the constraint submanifold M and we 
assume that there exist enough infinitesimal variations of this kind (that is, we are 
studying the so-called normal solutions of the vakonomic problem) . If M is locally 
given by the equations y A — ^ A (x l 1 y a ) = 0, for A = 1, . . . , m, we deduce that the 
infinitesimal variations must satisfy 

( lr I)^y A -$V,y a ))=0, X(x)=0, X(y) = 0. 

Note that if a G Adm([t a ,ti},M)% then 

(i v oX^ A (y A -^ A (x i ,y a ))o a = 

if and only if 

(ft ~ PaX OX* + dt 9y a + %a (313) 

+(^ + C>^. 
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Thus, if we consider our infinitesimal variations then 



d_ 

ds | s =o 



dL d^ J 



L ftl, QT 

= lie-;*; , dyA Qxi 



■X' 



dL 



-x c n 



dL d^ J 



It, 



\dx 



■p\x J 



Q y a a g yA g ya a 

dL 



X c n \dt 



BT 

—p l a X a + ^X c a ]dt 
dx l dy a a 



Now, let y A be the solution of the differential equations 



VA 



\dx l 



Vb 



dm 1 

dx l 



■)pA~yf(ci + ^ B ci B + y a ci a ), 



where 



dL 



d^ 1 



Qya A Qya 



(3.14) 



Then, from (|3.13p . we have that 

d 
dt 



d , - a s -a -a dL ■ - a , - „ d n< S> A 

-(y A X A ) = y A X A + y A X A = q~p\X a + y A p\X a — 



VA- 



dX a d^ A 
dt dy a 



-y A (C« + * B C« fl + y t>C« b )X~<^ + y A (CC + V B C A B 

+y c C A )X b - y b {C% + * B C b AB + y a C b Aa )X A . 
Using this equality, we deduce that 
d 

\dt Ki)/i ' dx' ,a dy a ' 



ds\ s =o J to 



to 



- d^ A dX a d^ A n i - d^ A 

-vaAX--^ - + y A (c« + **c« B + y b c« b) x^— 

+y b (C b M + * B C b AB + y a C b Aa )X A y A (C A + * B C A B + y c C A )X b )dt. 
Finally, using (|3.14[) and the fact that 



X c = 



dX a 
~df 



(C^ + ^ A C^ A + y b C^ b )X' 



we obtain that 

ds\s=oj to 



L{& s (t))dt = 



t dL dV J 
\dx l A dx l 



Pa 



-y l{ c: +* B c: B +y b cz b ) 



d / dL 
dt\dy a 

X a dt. 



d^ 



VA 



dy 



-) 



Since the variations X a are free, we conclude that the equations are 

pi + V A p A +y a P i a , 

Q^B, 



VA 



i dL 
\dx i 

d ( dl 
dl \ Thf 



VB 



VA 



dx'' 
q^a\ 

dy* . 



-)pA-y,(cl Q + ^ B ci B + y a ci a ), 



~ \dx l 



VA' 



d^ / 
dx % 



Pa 



- yi (c2 + * B c: B + y b c2 b ) 



dl d^ A 

for all 1 < i < in, 1 < A < m and m + 1 < a < n, with y a = — VA~^ — , that 

dy a dy a 

is, the vakonomic equations for the system (L, M) on t a : A — > Q (see (|3.6[1 ). 
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4. Examples 

4.1. Skinner-Rusk formalism on Lie affgebroids. Consider on a Lie affgebroid 
t/l : A — > Q a vakonomic system (L,M) with M = A, that is, a free system. In 
this case, Wq = A + ©q A and the Pontryagin Hamiltonian Hw ■ -A + ©q -A — > R 
is defined by Hw (<p,&) — ¥>(a) — L(a). Moreover, the precosymplectic structure 
(SlwoiV) on T A Wo is given by 

n Wo = (7pT 1 ,Wi)*^j l + d 7Aw °H Wo Arj and <q = (p~r 1; i/J'l^, 

where pr x : .A+ ©q .A — > j4 + is the canonical projection on the first factor, Tpi^ : 
T A (A+ ffi Q A) -> 7 A A+ is its prolongation and pr : : T A {A+ ffi Q A) -> A is the 
restriction of the projection A x T(A + @q A) — > .A (on the first factor) to the 
prolongation T^-A 4 " ©q A). In local coordinates, we have that 

Hw (x l ,yo,y a ,y a ) = yo + y a y a - L{x\y a ), 

and 

Then, the submanifold W[ C ,A + ©q .A is locally characterized by 

dL 

y = L(x\y a ) - y a y a , y a - — = 
and the vakonomic equations reduce to 

f ±* = ph + y a P^ 

\ Jtidr) = pla d^~ {c ^ + c ^ y(3) dy^' 

Thus, if (pr 2 )|vKi ■ — * -A is the restriction to W\ of the canonical projection 
on the second factor and 71 : / — > W\ is a solution of the vakonomic equations, 
then (pr 2 )|Wi °7i is a solution of the Euler-Lagrange equations for L. 

Note that in the standard case, that is, HA — J 1 t 1 this procedure is the Skinner- 
Rusk formulation for time-dependent mechanics (see [HE]). 

4.2. The 1-jet bundle of local sections of a fibration. Let r : Q — > R be 

a fibration and r^o : J 1 t — > Q be the associated Lie affgebroid modelled on the 
vector bundle 7r = (ttq)\vt '■ Vt —* Q (see Section HO]) . If (t,q z ) are local fibred 
coordinates on Q then {37, ^} is a local basis of sections of ttq : TQ — > Q. 

Denote by (i, (f, i, g 1 ) the corresponding local coordinates on TQ. Then, the (local) 
structure functions of TQ with respect to this local trivialization are given by 

C% = and p) = Sa , for i,j, k G {0, 1, . . . , n}. (4.1) 

Moreover, (i,^,^ 1 ) are the corresponding local coordinates on J 1 t. 

Now, let M C JV be a constraint submanifold such that t^oim : 3VC — » Q is a 
surjective submersion and L : J t — ► R be a Lagrangian function. Suppose that 
the constraint submanifold M is locally defined by the equations q A = \& A (t, q l , q a ), 
where we use the following notation (t, q l , q l ) = (t, q z , q A , q a ). 
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Thus, if we apply the results of the Section [3~T1 to this particular case, we recover 
some results obtained in [Tj. In particular, using (13. 6p and (|4.1|) . it follows that the 
vakonomic equations reduce to 



dL d^ B 

(hp 



PA = ~ T -VB p. A 



d/dL d^ A \ dL d^ A 



dt\dq a PA dq a ) dq a PA dq a ' 
q A = V A (t,q\q a ), 

where (t,q l ,p,pi) are the local coordinates on T*Q induced by the local basis {dt, 
dq 1 }. If these equations are written using the Lagrange multipliers (see Q3.7p ) then 
they coincide with the equations obtained in [33] . 

On the other hand, if (L,JA) is a regular vakonomic system on t\jo : J 1 r — ► Q, 
then the AV-bundle /ii : W\ — > W{ is locally defined by pi(t,q l 7 p,pi) — (t,q l ,pi). 
Thus, if h[,h'{ : W{ -> Wi are two sections of (inWi^y W{, 

h[(t,q\p t ) = (t,q % ,-H[(t,^,pj),pi) and h"(t,q\pi) = (t,q\-H"(t,q J ,Pj),Pi), 

then, from (|3.12p , we deduce that the vakonomic bracket {•, -} va k ■ T(mi) x T(mi) — * 
C°°{W[) is locally given by 

{h > ; h » }vak = d ^ dHjdHi 9H[dH'{ 

4.3. Optimal control systems as vakonomic systems on Lie affgebroids. 

Let ta : .A — > Q be a Lie affgebroid and C a fibred manifold over the state manifold 
7T : C — > Q. We also consider a section cr : C — > A along 7r and an index function 
I : C — ► R. The triple (Z, 7r, <t) is an optimal control system on the Lie affgebroid A. 

One important case happens when the section a : C — > .A along 7r is an embed- 
ding. In such a case, we have that the image M = cr(C) is a submanifold of A. 
Moreover, since a : C — > M is a diffcomorphism, we can define a Lagrangian func- 
tion L : M — >Rby£ = Z°er~ 1 . Therefore, it is equivalent to analyze the optimal 
control problem defined by (Z, 7T, a) (applying the Pontryagin maximum principle) 
that to study the vakonomic problem on the Lie affgebroid ta ■ A — > Q defined by 
(L,M). 

In the general case (when a : C — > A is not, in general, an embedding), we 
consider the subset d A C of the product manifold A x TC defined by 

3 A C = {(a, v) e A x TC | PA {&) = (Ttt)(«)}. 

Next, we will show that 3 A C admits a Lie affgebroid structure. Let r~ : 7 A C — > C 
be the prolongation of the bidual Lie algebroid .A of .A over the fibration 7r : C — > Q. 

defined by 



On the other hand, let <j) : 7 A C -> R be the section of (r|)* : (7 A C)* -> C 



</i(a,Xp) = lyi(a), for (a,X p )eT^C, with p E C. 

Note that if pr x : T^C — > A is the canonical projection on the first factor then 
(pr 1; tt) is a morphism between the Lie algebroids r~ : T^C — > C and : A —* Q 
and, moreover, we have that (pr 1 ,7r)*(l J q.) = 0. Since lyi is a 1-cocycle of : A — > 



T A ( 

the fact that (lyi),^ 7^ 0, for all x £ Q, we have that <A t .a C 7^ 0, for all p e C. 



Q, we deduce that <j> is a 1-cocycle of the Lie algebroid r~ : T C — > C and, using 
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In addition, it follows that 

= {(a,X p ) G TfC | U(a) = 1} = 3 A C. 

On the other hand, let Ty : 7 V C — > C be the prolongation of the Lie algebroid 
(V, [v]vj Pv) over the fibration 7r : C — > Q. Then, it is easy to prove that 
^{O} = T y O. Therefore, we conclude that 3 A C is an affine bundle over C with 
projection t\ : 3 A C — > C defined by rj(a, u) = 7rc(i>), where 7rc : TC — > C is 
the canonical projection. Moreover, the affine bundle : J^C — > C admits a Lie 
affgcbroid structure such that its bidual Lie algebroid is just (7 A C, [•, p ~) and 
it is modelled on the Lie algebroid Ty : 7 V C — > C (see Section HH]). 

Thus, we can consider the constraint submanifold M of the Lie affgebroid 3 A C 
defined by 

M= (J {(a,X p ) e#tf|<r(p) =a} 
pec 

and the Lagrangian function L : 3 A C —> R given by L = i °tJ. Then, (L, M) is i/ie 
vakonomic system associated with the optimal control system. 

If A = J x t is the 1-jet bundle of local sections of a fibration r : Q — ► ffi, it 
is easy to prove that the prolongation of J x r = TQ over 7r : C — > Q is just TC. 
Thus, a^C = {X G TC\dt(X) = 1} = J^t^tt), i being the usual coordinate on 
M. Under these identifications, the constraint submanifold is 

M={l£ TC | Ttt(X) = a(7rc(Jf ))}. 
Therefore, we recover the construction developed in PQ. 

Example 4.1. We consider the following mechanical problem (see [3ll4in5 l [2T | l29j). 
A (homogeneous) sphere of radius r > 0, mass m and inertia mk 2 about any 
axis rolls without sliding on a horizontal table which rotates with time-dependent 
angular velocity about a vertical axis through one of its points. Apart from the 
constant gravitational force, no other external forces are assumed to act on the 
sphere. The configuration space of the sphere is Q = M 3 x SO(3) and the Lagrangian 
of the system corresponds with the kinetic energy 

K(t,x,y;i,x,y,ui x ,0J y ,0J z ) = ^{mx 2 + my 2 +mk 2 (ujl + w 2 

where (u) x , u) y , u) z ) are the components of the angular velocity of the sphere. 

Since the ball is rolling without sliding on a rotating table then the system is 
subjected to the affine constraints: 

x — rujy = —Q(t)y, y + ruj x = £l(t)x, 

where f2(t) is the angular velocity of the table. Moreover, it is clear that Q = 
R 3 x 5*0(3) is the total space of a trivial principal 5'0(3)-bundle over R 3 and the 
bundle projection tt : Q — > R 3 is just the canonical projection on the first factor. 
Therefore, we may consider the corresponding Atiyah Lie algebroid TQ / SO(3) over 
R 3 (see [201 [22]). 

Since the Atiyah Lie algebroid TQ / 50(3) is isomorphic to the product manifold 
TR 3 x so(3) ^ TR 3 x R 3 , then a section of TQ/SO(3) ^ TR 3 x R 3 -» R 3 is a 
pair (X, u) , where X is a vector field on R 3 and u : R 3 — > R 3 is a smooth map. 
Therefore, a global basis of sections of TR 3 x R 3 — > R 3 is 
( d d d \ 

where u\, U2, M3 : R 3 — > R 3 are the constant maps Ui(t, x, y) = (1, 0, 0), U2(t, x, y) = 
(0, 1, 0) and u 3 (t,x,y) = (0,0,1)- 
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The anchor map Ptq/so(3) '■ TM 3 x R 3 — > TR 3 is the projection over the first 
factor and if [•, -Jtq/so{3) is the Lie bracket on the space T(ttq/sq^) then the 
only non-zero fundamental Lie brackets are 

[e3, e^TQ/soci) = es, [e4, e 5 ] T Q/so{3) = e 3, [es, e 3 ] T Q/so(3) = e 4- 

Moreover, : TR 3 x R 3 — > R given by <fr(t,x,y;t,x,y,u> x ,(J v ,u) z ) = t is a 1- 
cocycle in the corresponding Lie algebroid cohomology and, then, it induces a Lie 
affgcbroid structure over j4 = _1 {1} = R x TM 2 x M 3 . In addition, the affinc 
bundle t a : R X TR 2 x R 3 -> R 3 is modelled on the vector bundle t v : V — 
^{OJeIx TM 2 xl 3 ^! 3 and its bidual Lie algebroid A is just the Atiyah Lie 
algebroid TR 3 x R 3 . Note that the Lie affgebroid structure on A = R x TR 2 x R 3 
is a special type of Lie affgebroid structure called Atiyah Lie affgebroid structure 
(see Section 9.3.1 in [T7] for a general construction). Thus, (t, x, y; x, y, u> x , w y , oj z ) 
may be considered as local coordinates on A and V. 

It is clear that the Lagrangian function and the nonholonomic constraints are 
defined on the Atiyah Lie affgebroid A = R x TR 2 x R 3 (since the system is 50(3)- 
invariant). In fact, we have a nonholonomic system on the Atiyah Lie affgebroid 
A = R x TR 2 x R 3 (see [15] for more details). 

After some computations the equations of motion for this nonholonomic system 
may be written as follows 

x — rujy — —Q(t)y, y + ru> x = Cl(t)x, u> z — c, (4-2) 

where c is a constant, together with 

Now, we pass to an optimization problem. Assume full control over the motion of 
the center of the sphere and consider the cost function 

L(t, x, y; x, y, u x ,w y , uj z ) = - ((x) 2 + (yf) 

and the following optimal control problem: Given two points qo,qi G Q, find an 
optimal control curve (t, x(t), y(t)) on the reduced space that steer the system from 
qo and q\ and minimizes J Q | ((i) 2 + (y) 2 ) dt subject to the constraints defined by 
equations (14.2)) . 

Note that this problem is equivalent to the optimal control problem defined by 
the section u:R 3 xI 2 ^Mx TR 2 x R 3 along R x TR 2 x R 3 -> R 3 given by 

a(t,x,y;u l ,u 2 ) = (t, x, y; u 1 , u 2 , -{-u 2 + Q,(t)x), -{u 1 + fl(t)y),c) 

r r 

and the index function l(t, x, y; u x , u 2 ) = ^((u 1 ) 2 + (u 2 ) 2 ). Since a is obviously an 
embedding, we deduce the equivalence between both problems. 

A necessary condition for optimality of the problem is given by the corresponding 
vakonomic equations. In this case, we will denote by 

V 1 =x, y 2 = y, y 3 = w x , y A = uj v , y 5 = w z . 

Therefore, the vakonomic problem is given by the Lagrangian 



L(t,x,y iy \y 2 ,y 3 ,y\y 5 ) = ^((yr + (y 2 ) 2 ) 
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and the submanifold M denned by the constraints 



y 3 = 


* 3 (t, 


x,y,y\y 2 ) = ^-(-y 2 + fl(t)x), 


2/ 4 = 


* 4 (i, 


x,y,y\y 2 ) = -(y 1 + fi(t)y), 
r 


y 5 = 




^,2/,2/\?/ 2 ) = c. 


After some computations, we obtain that the vakonomic equations are 




m 




-^{y 1 + Q(t)y)y 5 + cy 4 , 




2/4 




-i(y 2 - Q(t)x)y 5 - cy 3 , 


< 


2/5 




^(y 1 + fl(t)y)y 3 - ^-(-y 2 + Q(t)x)y 4 , 




d ( 1 \ 

dt ( rv - ^ 




-n{t)y 3 , 




j t (ry 2 + y 3 ) 




-fl(t)y 4 , 




2/ 1 = ±, 




y 2 = y ■ 



Moreover, it is easy to prove that the vakonomic system is regular. Therefore, there 
exists a unique solution of the vakonomic equations on the submanifold W[ which 
is determined by the following conditions 



dL 



2/1 = JTT ~ VA-^-r = y - -2/4, 
oy L oy L r 

2 1 

y +-2/3, 



dL 
dy 2 

2/0 = L- y A ^ J 



dy 2 



l , 1 ^2 !/ 1 ^2 ^(*)/ s 

VaV = - o(2/i+ -2/4) - 77(2/2 2/3) {xy 3 + 2/2/4) -cy 5 . 

2 r 2 r r 



Thus, it follows that (t, x, y; yi,y 2 , 2/3, 2/4, 2/5) (respectively, (i, x, y; y , 2/1,2/2, 2/3, 2/4, 
2/5)) are local coordinates on W{ (respectively, on W\). Then, the local expression 
of the Hamiltonian Hy/ 1 is 



H Wl (t, x, y; 2/0, 2/1,2/2, 2/3, 2/4, 2/s) 



!/ 1 ^2 1, 1 ,2 

-2/0 - 7/(2/1 + -2/4 - 7; 2/2 - -2/3) 
2 r 2 r 



n(t) 



(^2/3 + 2/2/4) - cy 5 



and, in terms of the affine-linear part of the vakonomic bracket {-,-} V ak 

associated with the regular system (L,M), the vakonomic equations are 



2/i = {hi,yi}y l ak 

2/3 = {hl,y3}y ak 



n(t) 



1 



2/3, 2/2 = {/ll,2/2>Sife 



1 



fi(t) 



2/4, 



- 2/i + -2/4 + 2/5 + cy 4 



2/5 = e»i,Ws}„ ofc 



y 4 - {^i,2/4}£ifc = - ( -2/2 + -2/3 + ti(t)x ) 2/5 - cy 3 , 



2/i + -2/4 + fi(i)2/ 2/3 +( 2/2 2/3 - H(i)a; I y 4 



x = {hi,x}f ak = 2/1 + ^2/4, 2/ = {/H, y}vak = 2/2 - ^2/3- 



A 
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4.3.1. Optimal control of affine control systems. Let tji : A — ► Q be a Lie affgc- 
broid. Suppose that the constraint submanifold M of the vakonomic system is an 
affine sub-bundle 23 of A, that is, we have an affine bundle 23 over Q with associated 
vector bundle tu% '■ U-b — > Q and the corresponding inclusions is : 23 — > A and 
*[/ B : t^s ~~ ► V. Choose now a coordinate system adapted to this affine subbundle 
23. That is, take local coordinates (x l ) on an open subset U of Q and an local 
basis of r(ry), {e^, e a }, adapted to the decomposition V = U% U$, where U-b is 
an arbitrary complementary subspace. Thus, {eo, &A-, e a } is an local basis of r(rj) 
adapted to the affine subbundle 23, where lyi(eo) = 1. Denote by (a;*, y°, y A , y a ) 
the corresponding local coordinates on A and by {x % ,yo,yA,ya) the dual local co- 
ordinates on A + . Note that the equations which define 23 as an affine subbundle of 
A are y A = 0. 

The affine control problem given by the drift section eo and the input sections 
e a is defined by the following equation on Q, x l = p l + y a p l a , where the coordinates 
y a are playing the role of the set of admissible controls. 

Now, consider a function L : 23 — ► 1 as a performance index. The equations of 
motion of the optimal control problem defined by (L, 23) are precisely the vakonomic 
equations. In the selected coordinate system are: 

dL 

< VA = PAQ^i -y-y( C A0 + y aC Aa)> 

for all 1 < i < m, 1 < ^ < n, 1 < A < m and fn + 1 < a < n, with y a = 

Example 4.2. Consider a particle of unit mass in a planar inverse-square law 
gravitational field which has thrusters in the "x, y" directions (see [13]). Then, the 
equations of motion are: 

9l=«l, <Z2=w 2 , vi = -qi(qj + qj)~ 3/2 + v,!, v 2 = ~q 2 (qj + <?f)~ 3/2 + u 2 

defined on M = (R 2 — {(0, 0)}) x R 2 . The objective will be to drive the particle to 
a given circular orbit with minimum energy. Therefore, let us take L = \{u\ 

Now, choose a global basis of sections of T(K x M) — > R x M adapted to this 
affine control system: 

„ _ d 1 „ 9 1 „, d „ („2 1 _2\-3/2 d n (J2 , _2\-3/2 d 

„ d ~ d ~ d - d 

Ll — dqi ' C2 — dq 2 ' 3 ~ Stii ' C4 — 9i>2 

where {eo; 63,64} defines a affine subbundle of R x TM — > M. x M determining 
the initial affine control system. 

Applying the result developed in Subsection 14.3.11 and denoting by 2/3 = u% and 
2/4 = U2 then the equations of motion are now: 

qi=Vi, q_2 = v 2 
"1 = -5l(9l + <l 2 y 3/2 + Ui, v 2 = -q 2 {q\ + q 2 y 3/2 + u 2 
y x = - ( Ul {2q 2 - q 2 ) + 3u 2 q iq2 ) (q 2 + q 2 )- 5 ^ 2 
V2 = ~ (3u 3 q iq2 + Ui {2ql - q 2 )) (q 2 + q 2 )^' 2 

u\ = -y\ u 2 = -V2- 



A 
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5. Conclusions and future work 

Variational constrained Mechanics is discussed in the Lie affgebroid setting. We 
obtain the vakonomic equations and the vakonomic bracket associated with a con- 
strained mechanical system on a Lie affgebroid. The variational character of the 
theory is analyzed. Vakonomic systems subjected to afhne constraints are of special 
interest. Other examples are also discussed. 

In this paper we only consider normal solutions of the vakonomic problems. It 
would be interesting to extend the results of the paper for abnormal solutions. 
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